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Abstract. - The electronic structure of the corundum-type transition-metal oxides V2O3 and 
Ti203 is studied by means of the augmented spherical wave method, based on density-functional 
theory and the local density approximation. Comparing the results for the vanadate and the 
titanate allows us to understand the peculiar shape of the metal 3d a\ g density of states, which 
is present in both compounds. The oi s states are subject to pronounced bonding-antibonding 
splitting due to metal-metal overlap along the c-axis of the corundum structure. However, the 
corresponding partial density of states is strongly asymmetric with considerably more weight 
on the high energy branch. We argue that this asymmetry is due to an unexpected broadening 
of the bonding a\ g states, which is caused by hybridization with the e g bands. In contrast, 
the antibonding ai g states display no such hybridization and form a sharp peak. Our results 
shed new light on the role of the ai g orbitals for the metal-insulator transitions of V2O3. In 
particular, due to ai g -e g hybridization, an interpretation in terms of molecular orbital singlet 
states on the metal-metal pairs along the c-axis is not an adequate description. 



Based on the electronic level scheme of Castellani et al. [1], V2O3 was studied extensively 
as the canonical Mott-Hubbard system. Due to octahedral coordination of the metal sites, 
the V 3d states are split into lower i 2ff and higher e g levels. Trigonal lattice symmetry leads 
to further separation of the former into a\ g and states. Covalent V-V bonding along the 
hexagonal c-axis of the corundum structure results in bonding and antibonding a\ g molecular 
orbitals, which sort of bracket the levels. Electronic structure calculations confirmed the 
gross features of this scheme with t 2g states close to the Fermi level [2]. As proposed by 
Castellani et al., the bonding a\ g states are fully occupied, whereas the antibonding states 
shift to higher energies. This leaves one electron per V atom in the twofold degenerate e g 
orbital and leads to an S = 1/2 state, suggesting to use the one-band Hubbard model at 
half filling for describing the metal-insulator transitions (MITs) of V2O3. The stoichiometric 
compound undergoes an MIT at 168 K, leading from a paramagnetic metallic (PM) to an 
antifcrromagnetic insulating (AFI) phase; on doping with Cr or Al a paramagnetic insulating 
(PI) phase is found [3,4]. 

The model by Castellani et al. has been called into question by polarized x-ray absorption 
spectroscopy experiments, which point at an S = 1 spin state of the metal atoms [5]. The 
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lowest excited states miss pure e g symmetry but reveal remarkable a± g admixtures calling for 
a treatment beyond the one-band Hubbard model. LDA+U calculations explain the peculiar 
antifcrromagnetic order at low temperatures and point at an S — 1 configuration [6], which 
also results from the model calculations by Mila et al. [7]. Starting with the assumption of 
strong covalent bonding in the V-V pairs parallel to the c-axis, these authors suppose the 
intcrsite a\ g hopping matrix element to dominate. However, a recent study of the hopping 
processes in V2O3 found relevant matrix elements also between second, third, and fourth 
nearest neighbours [8]. LDA band structures show only minor response to the structural 
modifications occuring at the phase transitions of V2O3, whereas a correct description of 
the PM-PI transition has been obtained by a combination of LDA calculations with the 
dynamical mean field theory (DMFT) [9,10]. Nonetheless, confirming an early suggestion 
by Dernier [11] we have recently demonstrated that considering the rearrangements of the 
metal-metal overlap perpendicular to the hexagonal c-axis is important for understanding the 
MITs of V 2 3 [12-14]. 

Ti 2 03 is isostructural to corundum V2O3 and undergoes a gradual MIT without symmetry 
breaking lattice distortion between 400 K and 600 K. In contrast to V2O3, its low temperature 
insulating state is non-magnetic. As in V2O3, the a\ g states mediate strong bonding between 
Ti-Ti pairs in face-sharing octahedra parallel to the hexagonal c-axis, leading to bonding and 
antibonding a\ g bands bracketing the e g states. An insulating energy gap is expected between 
the bonding a\ g and e g states [15]. According to this picture increase of the c/a ratio of the 
corundum lattice constants with increasing temperature reduces the a\ g band splitting and 
promotes a collapse of the energy gap. Band structure calculations by Mattheiss [16] revealed 
a partially filled ti g complex of overlapping a\ g and e g bands at the Fermi level. Decreasing 
the c/a ratio reduces but does not eliminate the ai g -e g band overlap. To open the gap, an 
unphysically small Ti-Ti distance of 2.2 A parallel to the c-axis is required, which precludes 
a simple band explanation of the MIT. Recent cluster LDA+DMFT calculations assuming 
moderate Coulomb interactions among the t2 g orbitals reproduced the insulating state [17]. 

In this Letter we report on electronic structure calculations for both V 2 03 and Ti 2 03 
using the respective room-temperature crystal structure data. Our calculations i) result in 
a new interpretation of the bonding and antibonding a\ g states and ii) reveal striking dif- 
ferences in the hybridizations of the V 3d e g bands with the bonding and antibonding a\ g 
states, respectively. While hybridization of the e g bands with the bonding a\ g states is rather 
strong, leading to significant broadening of these bands, the antibonding states are of pure 
a\ g character. 

The present LDA band structure calculations are based on the scalar-relativistic augmented 
spherical wave (ASW) method [18,19]. Crystallographic data reported by Dernier for V2O3 
[11] and by Rice and Robinson for Ti2C>3 [20] were used. To model the crystal potential in 
the large voids of the open crystal structures, additional augmentation spheres were included. 
Optimal augmentation sphere positions and radii of all spheres were automatically generated 
by the sphere geometry optimization (SGO) algorithm [21]. The basis sets comprised V/Ti 
4s, 4p, 3d, and O 2s, 2p orbitals as well as states of the additional augmentation spheres. 
Brillouin zone integrations were performed using up to 2480 k-points within the irreducible 
wedge. 

As is typical for transition- metal chalcogenides with octahedral coordination, V2O3 shows 
three groups of bands in the vicinity of the Fermi level, extending from — 8.8eV to — 3.9eV, 
from — 1.1 eV to 1.5 eV, and from 1.9 eV to 3.8 eV, see Fig. ^ They originate from O 2p, V 
3d t 2g , and V 3d e g states, respectively. The right hand side of Fig. Q] shows the V 3d density 
of states (DOS) for the latter two groups of bands separated into its symmetry components 
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Fig. 1 - V2O3: weighted electronic bands and partial V 3d DOS. The length of the bars in the band 
structure indicates V 3d ai s contributions. 



ai g , e^, and e a g . V-0 hybridization results in additional O 2p admixtures in the energy range 
shown, which are stronger for the cx-bonding e a g states. Octahedral coordination of the metal 
atoms leads to almost perfect energetical separation of the tig and e g groups of bands. Small 
contributions of ti g states in the e g energy range above 1.9 cV are mainly due to the V-V 
anti-dimcrization along the c-axis, which shifts the metal atoms off the centers of their oxygen 
octahedra and results in a metal-metal distance much larger than the ideal value of 1 /6 of the 
c lattice constant. 

The band structure of the tig and e a g groups of bands is displayed on the left hand side 
of Fig. ^ where the length of the bars added to the bands is proportional to the respective 
a\ g contribution. Hence, bands without bars arc of and e a g character, respectively, in 
the lower and upper group shown. The band structure refers to the non-primitive hexagonal 
representation of the trigonal unit cell. While the dispersion of the e g bands is rather isotropic, 
the aig bands have an increased dispersion along the line T-A. Perpendicular to this line, i.e. 
along the paths T-M-K-r and A-L-H-A, the a\ g states are found mainly at the lower and 
upper boarder of the band complex at the Fermi level. This holds especially for the bands 
above leV in the T-M-K-r plane, which are well separated from the lower lying bands. 

The pronounced one-dimensionality of the a\ g bands leads to the characteristic shape of 
the ai g partial DOS. It suggests to interprete the peak at — 0.9eV and the double-peak at 
0.8/1.0 eV, respectively, as the bonding and antibonding states resulting from the V-V overlap 
across the octahedral faces. However, the weights of these peaks are far from equal but have 
a ratio of about 1:3, contradicting expectations based on a molecular orbital point of view. 
Due to the reduced weight of the lower a\ g peak the occupation of the twofold degenerate e g 
states is larger than one, explaining the experimental findings of Park et al. [5]. 

In order to understand these puzzling results we turn to titanium scsquioxidc, which has 
the same octahedral coordination of the transition-metal sites with octahedra neighbouring 
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Fig. 2 - T12O3: weighted electronic bands and partial Ti 3d DOS. The length of the bars in the band 
structure indicates Ti 3d a\ g contributions. 



along the c-axis interlinked via faces. The band structure and the partial densities of states 
of Ti2C>3 are displayed in Fig. [3 in the representation already used for V2O3 in Fig. From 
the structural similarity of the two compounds we expect the same gross features of their 
electronic structures. Indeed, we obtain for Ti2C>3 three groups of bands in the vicinity of the 
Fermi level, extending from -9.0 cV to -4.5 cV (O 2p), -0.9 eV to 2.3 cV (Ti 3d t 2g ), and 
2.4eV to 4.1 eV (Ti 3d e g ). Additional Ti 4s bands are found at energies higher than 3.9eV. 
These results agree well with the findings of Mattheiss [16]. 

The symmetry components of the Ti 3d DOS (right hand side of Fig. |2J) confirm the 
energetical separation of the t 2g and e g states as well as the t 2g admixture to the bands above 
2.4 eV due to the metal-metal anti-dimcrization along the c-axis. In the band structure shown 
on the left hand side of Fig. |2Jwe distinguish the isotropically dispersing e g bands from the 
aig states, which are highlighted by the bars attached to each band and display a rather 
one-dimensional behaviour. The latter causes the characteristic shape of the a\ g partial DOS 
with the pronounced peaks at — 0.7eV and 1.7 eV, which as before might be attributed to 
bonding-antibonding splitting due to metal-metal overlap along the c-axis. Yet, we are again 
faced with a decreased weight of the low-energy peak leading to the finite e g occupation. 

While the results for V2O3 and Ti203 are quite similar in general, we observe distinct 
differences at a closer look. They might give a first clue to the problem of the unequal weights 
of the dig peaks and pave the way for a better understanding of the role of these electronic 
states. To be specific, we point to the small gap at 1.2 cV in the partial DOS of the titanate, 
separating the peak at 1.7 eV from the lower lying bands. Such a gap is not observed for the 
vanadate. Since the a\ g partial densities of states both below and above 1.2 eV integrate to 
about one electron each, it is tempting to interpret the bands above this gap as antibonding 
and the whole of the a\ g bands between — leV and 1.2 eV as bonding. If this point of view 
were correct, the bonding a\ g states would thus extend over a very broad energy interval of 
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Fig. 3 - Hypothetical V2O3 (V-V bond lengths along Ch cx compressed to 2.51 A): weighted electronic 
bands and partial V 3d DOS. The length of the bars in the band structure indicates V 3d oi g 
contributions. 



more than 2 eV. 

This quite unusual situation should be related to the second observation growing out of 
the comparison of Figs.n an dEl namely the conspicuous downshift of the upper edge of the 

states. While in the vanadate these bands share edges with the whole of the a\ g states, 
there exists only a tiny contribution at energies above 1.2 eV in the titanatc. In contrast, the 

partial DOS displays a very sharp cutoff at this energy. Obviously, the e£ states are found 
only in the energy region of the bonding but not of the antibonding a\ g states. Going one 
step further, we infer from the similar shape of both partial densities of states below 1.2 eV, 
in particular from the common peaks near — 0.7eV, 0.2eV, 0.5eV, and l.OcV, a considerable 
amount of aig-e^ mixing. We thus attribute the large band width of the low-energy a\ g states 
to the strong ai g -e^ hybridization. Since the broadening is of similar size as the bonding- 
antibonding splitting of the a\ g band, we end up with the complex situation observed for the 
vanadate. 

In order to "prove" the just sketched scenario, we turn in the last step to another set of 
calculations for V2O3 with a hypothetical crystal structure. In this structure the V-V bond 
length along the c-axis has been reduced from 2.70 A to 2.51 A, while the symmetry of the 
corundum structure, the lattice constants, and the oxygen positions were retained. As a result 
of these changes, hypothetical V2O3 is characterized by a short intrapair V-V distance, hence, 
by increased intrapair bonding. Eventually, this will cause a stronger bonding-antibonding 
splitting of the a\ g states. As a consequence, we expect a much clearer separation of these 
states as compared to the results given in Fig. 2] The band shifts occurring on going from the 
real to the artificial structure of V2O3 will thus allow us to make a clear distinction between 
bonding and antibonding a\ g states for this material. 

The calculated results for the hypothetical V2O3 structure are displayed in Fig. using 
the same representation of the band structure and V 3d DOS as in Figs. Q and[21 While the 
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V 3d e g DOS of artificial V2O3 still resembles the results for real V2O3, we observe indeed 
distinct modifications especially for the a\ g DOS. The peak at — 0.8eV in Fig-^Jshifts to lower 
energies, and the peak at 1.0 eV (having a pronounced shoulder at 0.7 eV, see Fig. ^) splits 
into two peaks at 0.8 eV and 1.3 eV leaving a gap at about 1.0 eV. At the same time the upper 
edge of the e g states experiences a strong downshift and displays a sharp cutoff. Eventually, 
the band structures and densities of states of artificial V2O3 become more similar to those 
of Ti203. From the differences between real and artificial V2O3 we thus associate the a\ g 
peak around 1.3 eV with the antibonding states due to intrapair bonding. At lower energies, 
strong ai g -e g hybridization leads to a considerable broadening of the bonding a\ g branch, 
which extends from — 1.2 eV to 1.0 eV. Note that both the bonding and antibonding states 
integrate to equal weights. Transferring these findings to real V2O3, we conclude that the 
aig peak at 1.0 eV in Fig. ^ represents the antibonding states, whereas the region of bonding 
states extends from — 1.1 eV to « 0.9 eV and even includes the shoulder at 0.7 eV. 

In general, our findings agree well with the results of the recent LDA+U calculations 
by Ezhov et al. as well as Elfimov et al. for antifcrromagnctic and assumed ferromagnetic 
V2O3 [6,8]. The main difference between both approaches concerns the relative downshift of 
the e g states, which eventually opens the insulating gap observed for the AFI phase. However, 
the LDA+U calculations likewise result in a finite a\ g -e g hybridization for the bonding a\ g 
states whereas the e g admixture to the antibonding a% g states vanishes. Hence, while the 
LDA+U treatment leads to a somewhat changed scenario the basic mechanisms are still well 
described by the LDA calculations. 

Still, it remains an open question why the bonding a\ g states are subject to hybridization 
with the e g states, while the antibonding states retain their pure band character. Of course, it 
is tempting to relate this finding to the aforementioned metal-metal anti-dimcrization, which 
shifts the metal atoms along the c-axis off the center of their respective oxygen octahedra. 
Since the displacement is parallel to the principal axis of the a\ g orbitals, the resulting increase 
in d-p overlap mainly affects these states and is much less pronounced for the e g bands. The 
a\g states are thus found at slightly elevated energy, the difference of the centers of gravity 
amounting to « 0.3 eV for V2O3 [10]. As a consequence, the lower bonding part of the a\ g 
bands experiences more overlap with the e g states than the high-energy antibonding a\ g states. 

Nevertheless, a more convincing argument is based on the symmetry of the electronic 
orbitals involved, as well as the much different degree of overlap of the t2 g orbitals within 
the c-axis pairs. As a matter of fact, the e g orbitals arise as linear combinations of all five 
d states, except for the d^ z i- r i states which are identical to the a\ g states. In addition, the 
overlap of the orbitals within the vanadium pairs along the c-axis is close to negligible. 
As a consequence, these states are symmetric with respect to reflection about the midplane 
between the two atoms. This is different for the a\ g states, which do overlap along the c-axis 
this leading to the bonding and antibonding combinations. Since the bonding and antibonding 
a\g states are even and odd functions with respect to reflections, only the former may hybridize 
with the orbitals, whereas the latter do not overlap with these states. 

In conclusion, the electronic structure of the corundum- type sesquioxides V2O3 and Ti2 03 
is strongly influenced by a complex interplay of i) the formation of bonding and antibonding 
states arising from overlap of the metal 3d a\ g states across octahedral faces and ii) strong 
hybridization of the bonding but not the antibonding a\ g bands with the e g states. The 
hybridization disturbs the formation of molecular orbital singlet states by the a\ g orbitals of 
the vertical V-V pairs. Therefore, the occupation of the e g bands is increased this explaining 
the recently observed S = 1 configuration. As a consequence, failure of the model by Castellani 
et al. is mainly due to hybridization effects and not purely a consequence of the energy lowering 
of the e g states. In addition, since for symmetry reasons the ai g -e g hybridization affects only 
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the bonding a\ g states, the corresponding partial DOS assumes a strongly asymmetric shape. 
The observed strong 3D-coupling of the a\ g orbitals appears to be a primal feature of the 
electronic structure of these sesquioxides and therefore is of special importance for models 
dealing with the MITs in this compound. The findings for V2O3 strongly differ from those for 
VO2 and Ti407, where hybridization between the d\\ states, which are responsible for metal- 
metal overlap, and the e g states is suppressed and thus allows for a Peierls-type mechanism 
of the MIT [22,23]. 

* * * 
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